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DEFINITION OF SYMBOLS 


Symbol 

m 


P5 

V 


V > Pn 


m 


F 


s 


P6 

X 


t 

Pt 

IL 


d 

C. 


Definition 

mass of a meteoroid which at mean density, mean velocity, 
and mean impact angle will just puncture a given sheet 

logarithm of the puncture flux enhancement factor 

meteoroid velocity in kilometers per second 

meteoroid material density in grams per cubic centimeter 

log v and log p represent the population means of log v 
and log p p 

meteoroid mass in grams 

flux of incident meteoroids with mass equal to or greater 
than m 

value of log F s when m is one gram 

slope of log F s as a linear junction of log m 

angle of impact in radians with respect to the normal for 
a puncturing meteoroid 

effective thickness of a distributed puncture transducer 
in centimeters 

target subscript 

target material density in grams per cubic centimeter 

target hardness in Brinell units 

crater depth in centimeters 

projectile diameter in centimeters 

bulk velocity of sound in kilometers per second 

10“ 6 x Young's modulus, kilograms per square centimeter 

ductility, percent elongation in 2-inch gauge length 
at fracture 
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Symbol 

v 


k 2 

k 3 


k l 

0 

Pn 

\ 

0 


r 



i 


Pe 

f 
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DEFINITION OF SYMBOLS (Continued) 

Definition 

Poisson’s ratio 

a known function of target material parameters 

a parameter depending on material and dynamic properties 
of the projectile or meteoroid, thought to be 3 for 
dense projectiles 

an unknown function of meteoroid material parameters 
meteoroid puncture flux 

celestial latitude of radiant of surface normal 

celestial longitude of radiant of surface normal 

celestial latitude of earth's center seen from the 
satellite 

celestial longitude of earth's center seen from the 
satellite 

angular radius of the earth (meteoroid cross section) as 
seen from the satellite in radians 

.geocentric radial distance of the satellite in kilometers 
radius of the earth in kilometers 

inclination of meteoroid orbital plane to the ecliptic 
plane 

celestial latitude of meteor radiant 

cosmic weighting function for a photographic meteor 

meteor height above sea level at the point of maximum 
brilliance 

Opik's probability that a meteoroid in a given orbit 
will encounter the earth during one revolution of the 
particle 
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Symbol 


f 


f 


Z 

s 


x 2 


m x 


P 

mx 


R l 


R 2 

D 

7 

M 


Pr 


a 


z 
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DEFINITION OF SYMBOLS (Continued) 
Definition 


reciprocal of the apparent fraction of the circle of 
celestial latitude through a meteor radiant 

meteor zenith angle in radians 

constant scale factor to restore equality between the 
sarnie size and the sample sum of the weighting functions 
for photographic meteors 

angle of impact in radians with respect to the normal for 
a meteoroid impacting onto a surface 

mass of a meteoroid which at mean density, mean velocity, 
and angle of impact x will just puncture a given sheet of 
material 

probability that the mass m of a meteoroid is sufficient 
for puncture when the angle of impact is x 

the fraction of the potentially puncturing meteoroids 
from the exposed hemisphere which are intercepted by the 
apparent segment of the earth 

the apparent fraction of the earth 1 s disk 

angular separation of the center of the earth* s disk from 
the surface zenith in radians 

surface zenith angle (radians) of probable radiant 

angle at the center of the earth* s disk between the 
meridian of the exposed hemisphere and the celestial 
mer id ian 

celestial latitude of the probable radiant of the punctur- 
ing meteoroid 

celestial longitude of the probable radiant of the punctur- 
ing meteoroid 

half of the arc of the earth's disk (radians) subtended 
above the surface horizon 

angle (radians) between the centroid of the apparent disk 
and the center of the earth 

angle between celestial and surface meridians at the pro- 
bable radiant 
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THEORETICAL RELATIONSHIPS FOR METEOROID 
PUNCTURE EXPERIMENTS 

SUMMARY 


In unshielded meteoroid measurement satellite experiments where 
vehicle attitude information is available without information about the 
meteoroid closing velocity vector, the probable radiant and the puncture 
flux enhancement factor for sporadic meteors should be based on the 
assumption that velocity vectors for ambient meteoroids are isotropically 
distributed. Then, although half of the meteoroids which impact onto a 
flat transducer are more than 45 degrees from the normal, those within 
45 degrees puncture three times more effectively than the others, and 
half of the puncturing meteoroids are less than 33.5 degrees from the 
normal when the earth is below the surface horizon. This distribution 
also implies that the flux of puncturing meteoroids is five times greater 
than the flux of incident meteoroids with mass equal to or greater than 
the mass of the almost -puncturing meteoroid with the average values 
of density, velocity, and impact angle. 

The formulas for the celestial coordinates of the probable radiant 
of a puncturing meteoroid have been corrected for earth-shielding. This 
correction is important near the earth even under nominal circumstances; 
e.g., when the orbit height is 500 kilometers, the radius of the earth’s 
effective disk (as seen from the vehicle) is 70 degrees. In that case, 
when the elevations of the earth’s center are 0, 15, 30, 45, and 60 
degrees the probable radiants are deflected 5, 10, 20, 39, and 80 degrees, 
respectively, by the earth-shielding effect. 

I. INTRODUCTION 
Purposes 

The main purpose of this analysis is to show how attitude 
information for an appropriate meteoroid puncture counting satellite 
can be used to infer the probable radiant of a puncturing meteoroid even 
when the earth partially obscures the celestial hemisphere to which the 
puncture transducer surface is exposed. A secondary purpose prerequisite 
to the main purpose is the derivation of the probability density function 
for the angle of impact of a puncturing meteoroid. A biproduct of the 
secondary purpose is a more accurate determination of the puncture flux 
enhancement factor; i.e., if a particular sheet of material is just 
puncturable by a meteoroid of mean density, mean velocity, mean impact 
angle, and mass m, then the mean flux of puncturing meteoroids exceeds 
the mean flux of incident meteoroids with mass equal to or greater than 
m by (3 S order of magnitude [1]. The more accurate determination of (3 5 
is a secondary purpose of this analysis. 



Methods and Assumptions 

In most of this analysis, it is not necessary to know the true 
mean value of parameters such as the logarithms of the closing velocity 
and the material density p of the meteoroids, but only the variation with 
respect to the unknown true mean. All indicated logarithms are to base 
ten. The following assumptions seem appropriate: 

1. All puncturing meteoroids are considered to enter cis lunar 
space from outside the earth-moon system with appreciable velocity. 
Meteqroid geocentric velocity at any point in cislunar space is suf- 
ficiently higher than vehicle orbital velocity that it is essentially 
the closing velocity. The logarithm of velocity, v, when v is in kilom- 
eters per second, for the population of incident meteoroids, is essen- 
tially normally distributed with mean and standard deviation given by 
equations (1) and (2), respectively [1], 


log v = log 26.7 = 1.43, 


( 1 ) 


o\ = 0.18. 
log v 


( 2 ) 


2. Most of the puncturing meteoroids in all present experi- 
ments are dust balls instead of having the more substantial structure 
of meteorites. Although the mean density, p , for the population of 
dustball meteoroids is known only to within a factor two probable error, 
the logarithm of density, p p , where the density is in grams per cubic 
centimeter, for the population of incident meteoroids, can be considered 
to be normally distributed with mean and standard deviation given by 
equations (3) and (4), respectively, [1]. 


log p = log 0.44 ± 0.30 = - 0.35 ± 0.30, (3) 

P 

ex, = 0.44. (4) 

lo s Pp 


3. It is assumed that the puncture transducers are exposed to 
a celestial hemisphere, and that the probable radiant of the velocity 
vector of a puncturing meteoroid is determined with sufficient accuracy 
by neglecting any systematic variations in meteoroid physical or statis- 
tical parameters over that part of the apparent celestial hemisphere 
which is not subtended by the earth. Most of the punctures are expected 
to be from sporadic meteors rather than from stream meteoroids [1]; but 
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the known occurrence of a major meteor stream is an exception to the 
isotropic flux which is otherwise assumed. If a puncture is concurrent 
with a well-defined meteor stream with unobstructed radiant within the 
median angle of deviation from the surface normal radiant, then the 
probable radiant should be considered to be that of the concurrent 
meteor stream. 

4. At any particular distance from the surface of the earth, 
the mean flux, F s , of meteoroids per unit area per unit time with mass 
equal to or greater than m grams is considered to be an exponential 
function of mass m; i.e., [1], 


F = 10 Ps m^ 2 
s 


(5) 


where (3 6 is a constant. In a model for cis lunar meteoroid impact recently 
suggested as a .criterion for current purposes by Dalton [1], is given 
as a function of distance from the surface of the earth with range 
-1 . 34 ^ - 1 in cislunar space and with an effective value of -1.23 

for the orbit of the Pegasus -A meteoroid detection satellite. To deter- 
mine the distribution of the angle of impact, x, for the population of 
puncturing meteoroids, the numerical value of for the Pegasus-A mis- 
sion will be assumed throughout cislunar space; i.e.. 


® - 1 . 23 . ( 6 ) 

5. In Dalton’s [2] model for the puncture of a homogeneous 
metallic wall by hypervelocity impact of a spherical projectile, the 
relation between the minimum sufficient mass, m, of the projectile, 
its velocity, v, and the angle of impact, x, can be expressed as 


log m = 1.9 + log(p 3 p t H t ) - (3/2) (log v + log cos x) - log p , 

(7) 

where p is the density of the projectile and where p t , H t , and p are 
the deniity, hardness in Brinell units, and thickness in centimeters, 
respectively, for the punctured wall. Equation (7) is used in the 
present analysis only in the derivation of the distribution of the 
angle of impact, x, for the population of puncturing meteoroids. 
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6. Under contract number NAS7-219 to the Martin-Mar ietta 
Corporation for a survey and multivariate analysis of available labora- 
tory hypervelocity impact data, Reismann, Donahue and Burkitt [3] found 
that in the highest of three velocity regimes the multiple correlation 
coefficient is neatly maximized (to 0.989) by representing the relation 
between projectile and target parameters for impact at normal incidence 
as follows: 


P Q /d = 2.5(v/C t ) 






V.J 


where 


p Q = crater depth 
d = projectile diameter 

v = closing velocity in kilometers per second 
t = target subscript 

C t = bulk velocity of sound in kilometers per second 
E = 10“ 6 x Young's modulus, kilograms per square centimeter 
e = ductility, percent elongation in 2-inch gauge length at fracture 
v = Poisson's ratio. 

The exponents in equation (8) are the result of some rounding off of 
decimals in the published [3] values. The relation between the thick- 
ness, p, of a puncturable sheet and crater depth, p Q , in a thick target 
has been treated analytically by Andriankin and Stepanov [4] and by 
Herrmann and Jones [5] and stochastically by Dalton [2, 6]. The results 
are problematical; but for meteoroid impact, the following relation 
would seem to be the most appropriate assumption’: 

p = 1.59 p 0 . (9) 


Dalton [1] showed that when equations (1), (3) and (9) are assumed for 
meteoroids, the compatibility of equations (7) and (8) for aluminum 
2024-T3 implies that the mass, m, - of a meteoroid with mean log density 
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(log Pp = log 0.44), mean log velocity (log v = log 26.7), and mean 
impact angle (x 2 = jt/4) - which will just puncture a particular metallic 
wall is 


fn = 10 1O,y85 ('DC °.50 E 1.31 -u S -° p 0.43\3 

VF t t t t ' * 


( 10 ) 


It is assumed that equation (10) is superior to equation (7) for those 
purposes where it can be used; however, the exponent "3" and consequently 
also the exponent "10.985" in equation (10) are problematical for mete- 
oroid impact, as well as the corresponding values in equation (7), and 
should be treated as unknown constants, say k 3 and k x , respectively, 
when incident flux parameters (p 6 and in equation (5)) are to be 
inferred from puncture flux parameters (see Section II below). 

1 1. RELATION BETWEEN PARAMETERS FOR IMPACT AND 

PUNCTURE FLUXES 

Corresponding to the meteoroid impact flux, F g , in equation (5), 
there will be a puncture flux, 0 , per unit area per unit time, through 
a particular exposed wall or distributed puncture transducer, which can 
be represented by 


0 = io^s+Ps 


( 11 ) 


where (3 6 and are the same constants used in equation (5) to express 
the incident (or impact) flux, F s , the constant, p 5 , is the 
logarithm of the puncture flux enhancement factor to be derived in 
Section IV, and m is the nominally puncturable meteoroid mass which is 
now predicted by the improvised equation (10). 

By the suggestion at the end of Section I about the problematical 
exponents in equation (10), an expression for m alternative to equation 
(10) is 


m 


10 


k -L+k 2 k 3 



( 12 ) 


where and k 3 are unknown functions of the projectile material and 
dynamic characteristics, to be determined for meteoroids by satellite 
experiments, and where k 2 is a function of the target material param- 
eters as follows: 


(13) 


k 2 = 0.50 log C t + 1.31 log E t + 8.0 log Vfc + 0.43 log e fc . 
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By equations (11) and (12), one can write the following relation 
between log 0 and (k 2 + log p), both of which (by equation (13)) have 
known values for each of the puncture transducers of different material 
and effective thickness in the flight experiments: 


log 0 = (P6 + P 5 + k 10 2 ) + ( k 2 + log p). (14) 


Dalton's [1] calculations for the Explorer XVI satellite imply 
-5.417 and -5.697 for log 0 for the 1-mil and 2-mil puncture sensors, 
respectively, and -6.175 and -5.904 for (k 2 + log p) for the same 
respective sensors. The 2-point solution to equation (14) for the 
Explorer XVI orbit would be -1.20. 


log 0 = -11.80 - 1.033 (k 2 + log p). 


(15) 


The resulting slope in equation (15) would imply by equation (14) that, 
if 3 were the correct value of k 3 , then the value for j3 2 would be -0.344. 

The value of p 2 which Hawkins and Upton [8] found for photographic meteors 
is -1.34; and by Dalton's [1] model the effective value of p 2 for the 
Explorer XVI orbit would be -1.20. 

Extrapolation of the puncture results for Explorer XVI by equation 
(15) would indicate a higher puncture rate for thicker structures than 
one would want to anticipate for design purposes without verification 
by further puncture experiments. The results in equation (15) are 
not established with much confidence because the numbers of punctures 
(44 and 11 for the 1-mil and 2-mil sensors, respectively) are statis- 
tically small samples, especially when the extreme diversity of the 
population of incident meteoroids is considered. More pertinently, one 
should anticipate that the relation between log 0 and (k 2 + log p) in 
equation (14) may not be linear. 

A sufficiently accurate value for p 5 will be derived in Section IV; 
but this leaves more unknowns (4) in equation (14) than conditions (2) 
which can be applied toward their determination from Pegasus-type data. 

If values for (3 6 and p 2 can be assumed from the statistics and physical 
theory of meteors, then the meteoroid impact parameters k x and k 3 can 
be found as constants (if equation (14) is linear) or as functions of 
(k 2 + log p) (if equation (14) is nonlinear). But, because the impact 
parameters k 1 and k 3 are not known, it is not even possible to determine 
the value of 3 6 in equation (14) by assuming a value for p 2 , or vice 
versa. Therefore, by equations (5), (11), (12), and (14), it follows 

that the nominally minimum sufficient mass for puncturing meteoroids 
cannot be determined from meteoroid puncture data alone, unless the incident 
flux, cumulative with respect to mass, is known as a function of the mass. 



III. DISTRIBUTION OF ANGLE OF IMPACT FOR INCIDENT 
METEOROIDS ISOTROPICALLY FROM AN APPARENT 
HEMISPHERE ONTO A PLAIN OR CONVEX SURFACE 


Appropriate puncture experiments should give satellite attitude and 
punctured sensor identification at each meteoroid puncture event (see 
Reference 9 for a discussion of the Pegasus-A experiment) . This informa- 
tion can be used to determine the celestial latitude, f3 n , and celestial 
longitude, A n , of the radiant of the normal to the punctured surface, the 
corresponding celestial coordinates (f3 p , A p ) for the center of the earth 
(as seen from the punctured satellite), and the angular radius G of the 
earth's effective meteoroid encounter cross section seen from the 
satellite; i.e., approximately, 

0 = sin -1 [ (r e + 100)/r], (16) 

where r is the geocentric radial distance of the satellite in kilometers 
and r e is the radius of the earth in kilometers. Therefore, when the 
probability distribution of the radiants of incident meteoroids is assumed, 
and when the relation between the probability distributions of the angles 
of impact of incident meteoroids and of puncturing meteoroids is estab- 
lished (see Section V), the celestial coordinates (0m>An) of the probable 
radiant of the puncturing meteoroid can be determined as functions of 
the values of A n , p p , A p , and 0 at the time of puncture. In this 

manner, any large sample of punctures will give a distribution of pro- 
bable radiants, although the resulting distribution may be somewhat more 
widely dispersed than the unknown population of true radiants. 

The dispersion of the probable radiants of the puncturing meteoroids 
in appropriate experiments will have two spurious components because of 
some uncertainty in the satellite attitude information and because the 
meteoroid velocity vector in on-board coordinates is not measured but 
must be inferred probabilistically. Nevertheless, The resulting infor- 
mation is of considerable interest because the other three sources of 
such information are somewhat problematical and they involve different 
ranges of meteoroid mass than that for which puncture experiments are 
designed. One such source of information is the interpretation of 
isophotes of zodiacal light from small dust particles in interplanetary 
space. Other inferences can be attempted from the distribution of the 
orbits of comets and from the effect of the combined action of several 
physical processes on the debris from disintegrating comets. More 
intensive studies of the distribution of meteoroid radiants have been 
based on the visual, photographic, and radio meteor data; but several 
sources of observational bias complicate these analyses and may pre- 
judice the results. 
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Without specific contrary information, one could assume that the 
distribution of the closing velocity vectors is isotropic and therefore 
that the probability density for meteor radiants is uniformly distributed 
over the exposed celestial hemisphere. Then the probability density 
function f(p e ) for the celestial latitude p e of the radiants in the 
interval 0 ^ p e ^ jt/ 2 would be cos f3 e ; and, by integration, the mean 
(expected value) and median for (3 e would be (jt/ 2 - 1) and it/ 6 radians 
(32.7 and 30.0 degrees), respectively. 

Lovell's [10] results for the distribution of sporadic visual 
meteors, corrected for horizon zenith attraction and twilight effects, 
have been interpreted by Davison and Winslow [11] to indicate that 
approximately half of the sporadic meteoroids have orbital inclinations, 
i, within ±15 degrees from the ecliptic. (Orbit inclination, i, in the 
interval 0 ^ i ^ jt/2, and radiant celestial latitude, p e , are practically 
interchangable here.) But the masses of the meteoroids had to be suf- 
ficient for duplicate visual observations (see Lovell's [10] description 
of Prentice's technique for visual meteors), and therefore several orders 
of magnitude higher than the largest mass one can practically expect to 
detect in present experiments. Then, by Hawkins [12] models for the 
distributions of meteoroids of cometary origin (most of those encountered 
in present experiments) and of asteroidal origin, the 2-observer visual 
meteor sample might be expected to have angles of inclination, i, 
distributed somewhat more like the asteroids. Allen [13] gives the 
arithmetic mean for the inclination of the orbits of the asteroids as 
approximately 10 degrees. 


Assuming an isotropic distribution of the radiants of meteoroids 
(encountered by an observer in orbit about the sun) is not equivalent 
to assuming completely randomly distributed meteoroid orbits (with 
inclination, i, in the interval 0 ^ i ^ jt/2). Briggs [14] indicates 
that sin i would be the probability density function for the inclina- 
tion, i, of completely random orbits (with all meteoroids considered, 
i.e., not just those detected by a moving observer). Then, by integra- 
tion, the mean (expected value) of i would be one radian (57.3 degrees). 

A completely random distribution of orbits may be indicated for 
long period comets; but the debris from the disintegration of comets 
may have orbits which are distributed less randomly. Davison and 
Winslow [15] describe this situation as follows: "The zodiacal dust 

particles are thought to be more highly concentrated in the plane of 
the ecliptic. They have two components of drift as they orbit the sun, 
one toward the plane of the ecliptic and one toward the sun [16]. The 
drift toward the plane of the ecliptic results from the attraction of 
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the planets. The drift or spiraling in toward the sun is a result of ... 
the Poynting-Robertson effect (Lovell [10] and Wyatt and Whipple [17]). n 

By establishing isophotes of the zodiacal light and interpreting 
them in relation to the Poynting-Robertson effect and data from observed 
meteors, Briggs [14] found that the probability density function f(i) 
for meteoroid orbLt inclination, i, (0 ^ i ^ jt/2) is 


f(i) 


(1.66 - 0.66i) sin i, for aphelion > 6.6 astronomical units 
— 6 i 

37e~ sin i, for aphelion S 6,6 astronomical units. 


Then, by integration, the mean (expected value) of the inclination, i, 
of the orbit of any high-aphelion particles (whether or not encountered 
by a moving observer) would be 1.66 - 0.66 (jr - 2) = 0.907 radians 
(52 degrees). The corresponding mean i for low-aphelion particles 
would be 12/37 radians (18.6 degrees). In this analysis, Briggs indicated 
that particles in the 1 to 50 micron range play the dominant role in the 
zodiacal light; and Ehricke [18] says that, according to Takakubo, the 
radius of the particles most effective is about 20 microns. Also, by 
the relation which Briggs indicated between particle radius and material 
density, a spherical part'icle with a radius of 20 microns would have a 
mass of 10" 8 * 2 grams. 

The author has completed only the first part [26] of his planned 
four-part analysis of the 285 photographic meteors common to the two 
tables of data for random samples of sporadic meteors published in 1958 
and 1961 by Hawkins and Southworth [19, 20]. Some preliminary results 
were given in Reference 1; e.g., the sample range for meteoroid mass, m, 
in grams (above the atmosphere) is 10~ 2 * 2 ^ m ^ 10 1- °, with a mean for 
log m of -1.0. Figure 1 shows the cosmically weighted cumulative 
distribution of celestial latitude, p e , with an arithmetic median f3 e of 
25 degrees for the positive values of p e . Asterisk and circle points in 
Figure 1 represent positive and negative values of g e , respectively. 

Figure 1 was obtained by weighting each meteor datum by 


f = 8 h- 2 v-3 /a P" 1 f. e“°- 18Z (17) 

m r 

where 1^ is the meteor height above sea level at the point of maximum 
brilliance, v is the closing velocity, P is Opik's (see Whipple [21]) 
probability that a meteoroid in a given orbit will encounter the earth 
during one revolution of the particle; ff is the reciprocal of the 
apparent fraction of the circle of celestial latitude through the 
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meteor radiant, e is the base of natural logarithms, Z is the meteor 
zenith angle in radians, and s is a constant scale factor which restores 
equality between the sample size and the sample sum of the weighting 
functions. The constant, -0.18, in the zenith angle weighting factor 
in equation (17) was chosen to minimize the sum of the squares of the 
differences between the cumulative distributions for positive and 
negative p e below 42 degrees in Figure 1. (The zenith angle weighting 
factor is an approximate alternative to weighting with respect to 
meteor celestial longitude.) It is not yet known to what extent the 
results in Figure 1 may be changed by deleting Opik's earth encounter 
probability, P, from the weighting function, f (equation (17)); but 
the corresponding results would be more appropriate for the present 
consideration (observer subject to the earth 1 s orbital motion). 

In an analysis of the distribution of sporadic radio meteor radiants 
corrected for antenna selectivity, Elford, Hawkins, and Southworth [22] 
show contour diagrams of radiant density over the celestial sphere, for 
each of the first eight months of 1962. Also they show corresponding 
results for the entire period both before and after correction for 
velocity selection effects. The monthly mean results, not corrected 
for velocity selection effects, show much shifting of the contour diagrams 
of radiant density from month to month. The mean results for the entire 
period undergo major changes when corrected for velocity selection 
effects; e.g., the concentration of radiants around the apex of the 
earth's way disappears, and the previously supposed scarcity of radiants 
near the antipex becomes questionable. Correction for velocity selec- 
tion effects is important for the present purpose because the incident 
flux of radio meteors is proportional to the 4f3 2 power of the velocity 
[22]; whereas, puncture flux through a given structure is thought to 
be proportional to the 1.50 2 power of the velocity [1]. The parameter 
(3 2 is the slope of log F s as a function of log m, and is thought to 
have an effective value between -1.34 [8] and -1.0 [22] for meteor data. 

The main features of the preliminary results for the mean radiant 
distribution of sporadic radio meteors, corrected for velocity and other 
selection effects, are that radiants are more densely concentrated near 
the sun and antisun directions (about ±80 degrees latitude from the apex 
of the earth's way) and that radiant density decreases by a factor of 
about three within about 25 degrees from the center of each of the two 
main concentrations. Elford, Hawkins, and Southworth [22] indicate that 
these results are for meteoroids with mass not less than 1.3 x 10" 6 grams. 

The dispersion in the sample of probable radiants from puncture 
experiments will have a spurious component due to any error in the dis- 
tribution of radiants which must be assumed in determining each of the 
probable radiants. It would be prohibitively tedious to apply a dis- 
tribution with intricate variations with respect to celestial coordinates 
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and season of the year. For any distribution with respect only to 
celestial latitude, it is not apparent from the discussion so far in 
this section that any improvement over the isotropic distribution could 
be found. Therefore, for simplicity and because the radiant activity 
within ±60 degrees from the antisun is not yet established with cor- 
rections for velocity selection effects for radio meteors, any varia- 
tions in the distribution of radiants with respect to celestial longitude 
will be ignored also in the rest of this analysis. 

Particles incident isotropically from an exposed hemisphere onto 
an element of a plane surface give a sin 2 x 2 probability density function 
for the angle of impact, x 2 , with respect to the normal; i.e.. 


f(x 2 ) = sin 2x 2 . 


(18) 


This result (equation (18)) is evident intuitively because that fraction 
of the celestial hemisphere which contributes radiants within the angular 
increment, dx 2 , varies as sin x 2 , and the projection of the element of 
area onto the plane normal to any radiant direction varies as cos x 2 . 
Then, by integration, both the mean and the median for the angle of 
impact, x 2 , for all meteoroids incident isotropically onto a plane sur- 
face, are equal to jt/4 radians. 

IV. PUNCTURE FLUX ENHANCEMENT FACTOR 


Equation (7) gives the logarithm of meteoroid mass, m, as a linear 
function of the logarithms of impact velocity, v, and meteoroid density, 
p . But, in Section I, log v and log p are described as essentially 
normally distributed and statistically independent random variables. 

Then log m in equation (7) is a normally distributed random variable 
with mean and variance, by equations (1) through (4), as follows: 


log m x = - 1.9 + log(p 3 p t H t ) - (3/2) log v - log p p - (3/2) log cos x 

(19) 


o 


2 

log m 


(3/ 2) 2 a 


2 

log V 


0*7 

!og P p 


(0.52) 2 . 


( 20 ) 
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By taking the differential of equation (5), the number of meteoroids, 
-dF g , impacting per square meter per second with masses between m and 
m + dm is 


-dF s = -10^ 6 dm. 


( 21 ) 


Then, by equations (18) and (21), the flux of impacting meteoroids, 
within both an increment of mass dm and an increment of impact angle 
dx 2 , is 


-dF (s in 2xo) 
s ^ 


dx 2 = -10 Ps m P2_1 


sin 2x 2 dx 2 dm. 


( 22 ) 


By equation (22), the incremental number of punctures 60, per 
square meter per second, due to meteoroids in the increment of impact 
angle dx can be expressed by 


60 


r -10 Pb B 2 P m^ 2 " 1 
J mx 



sin 2x dx, 


o 


(23) 


where P^ is the probability that the mass, m, of a meteoroid, which 
is incident with an angle of impact, x, is sufficient for puncture 
(i.e., P mx is the cumulative distribution of the normally distributed 
variable log m as defined by equations (7), (19), and (20)). The sub- 
script has been dropped from x 2 in equation (23) to indicate that the 
population of puncturing meteoroids is now considered rather than the 
population of impacting meteoroids. 


The tangent line at the point P mx =1/2 for the curve relating 
P^ and log m is a very close approximation to the curve in the interval 
0.25 ^ P^ ^ 0.75, and its extension throughout the interval 0 ^ Pmx ^ 1 
is close enough for the present purpose; i.e., 


mx 


n r: . ( \ ^ 2jt 

= 0 for m < exp 1Q (log m^ - ~ 2 ~ 


log m 


= 1 for m > exp 10 [ log \ + ^ CT log m 


J~2n , 


/log m - log m 

- <2„)-l/2( 

log m 


+ otherwisej 


(24) 
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Therefore, by using the expressions from equation (24) in evaluating 
the ingegral in equation (23), it follows that the differential puncture 
flux, 60, corresponding to the differential impact angle, dx, is 


60 = 10 


P6 



2* p 2 


log m 


717 


10 


p 2 a 


logm 


- 10 


TST" 0 

2 “ 2 logrriN _(3 2 


A 


sin 2xdx. 


(25) 


In equation (19), by replacing log cos x with the mean of log cos x 2 , 
say log cos x 2 , one also replaces m x (the mass of the nominally puncturing 
meteoroid at impact angle x) with m (the mass of the nominally puncturing 
meteoroid in equation (11)). Then, by equation (18), 


*/2 

log cos x 2 = J' (log cos x 2 ) sin 2x 2 dx 2 = - (1/2) log e; (26) 
o 


and, by equations (19) and (26), 


log = log m - (3/4) log e - (3/2) cos x. 


(27) 


One finds, by substituting the antilogarithm of the right side of equa- 
tion (27) for m x in equation (25), and by using the numerical values of 
and o^ Q m from equations (6) and (20), respectively, that 


60 = 10 P6 m^ 2 (9.69 sin 2x) (cos x) 1 - 84 dx. (28) 


Then, by integrating the right side of equation (28) over the interval 
0 ^ x ^ it/ 2 , 


0 = (5.05) 10^ m^ 2 = A + °- 7 ° rr7 2 . (29) 
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Therefore, by equations (5) and (29), the puncture flux enhancement 
factor is five, instead of four [1], and g 5 in equation (11) is 0.70 
instead of 0.60. 


V. DISTRIBUTION OF ANGLE OF IMPACT FOR PUNCTURING 
METEOROIDS ISOTROPICALLY FROM AN APPARENT HEMISPHERE 
ONTO A PLAIN OR CONVEX SURFACE 

The probability density function f(x) of the angle of impact, x, 
for the population of meteoroids which actually puncture a given sheet 
of material exposed to an isotropic flux, is, by equations (28) and (29), 


f(x) = (60/ dx) / 0 = 1.92(cos x) 1 * 84 sin 2x. (30) 


Then, by integrating equation (30), the probability that a puncturing 
meteoroid had an angle of impact of not more than x radians from the 
normal is 1 - (cos x) 3,84 . Therefore, the median and third quartile of 
impact angle for puncturing meteoroids are 33.5 and 45.8 degrees, 
respectively, from the normal. 

VI. PROBABLE RADIANT FOR A PUNCTURING METEOROID FROM 
A HEMISPHERE PARTIALLY OBSCURED BY THE EARTH 


The fraction of the exposed celestial hemisphere which has the 
meteoroid radiants corresponding to a differential increment of impact 
angle, dx, varies as sin x. Then, by equation (30), the number of 
meteoroids which puncture a hemispherically exposed flat transducer and 
which have radiants within a constant differential solid angle in the 
vicinity of a point on a meridian of the exposed hemisphere is propor- 
tional to (cos x) 2,84 . 

The hemispherical mean radiant density for the puncturing meteoroids 
is not expected to differ appreciably from the radiant density correspond- 
ing to the median impact angle for puncturing meteoroids (for which 
(cos x) 2 * 84 equals 2“ 84 ' 3, 84 ) ; i.e., (1.67)” 1 ). 

The apparent disk of the earth with angular radius, 0, in equation 
(16) subtends a solid angle of 2jt(l - cos 0) steradians, (see Refer- 
ence 23). Then, of the total number of meteoroids which would puncture, 
the fraction R ± obscured by the apparent fraction R 2 of the earth* s disk 
is sufficiently approximated by 
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R x = 1.67 (cos x ) 2 * 84 (1 - cos 0 ) R 2 , (31) 

where x is the separation of the centroid of the apparent segment from 
the surface zenith. 

The separation D of the center of the earth's disk Op,^p) from 
the surface zenith (£3n*A n ) is related to the given coordinates through 
the formula for the haversine of D (see Reference 25): 

hav D = hav (8 - 6 ).+ cos 6 cos ft hav (A - A ). (32) 

VH p l rr K n K p n p y K J 

By definition, the haversine and cosine of any angle A are related by 
hav A = (1 - cos A)/2. (33) 


Then, by equations (32) and (33), D is the smallest positive value in 
radians satisfying 


cos D = cos (ft - ft ) - [1 - cos (A - A )] cos f3 cos p , (34) 

K p n n p n p 

where D is in the interval 0 ^ D ^ jt. 

When D ^ jr/2 + 0, the earth is below the surface horizon, both 
R 2 and R-j_ in equation (31) vanish, and the probable radiant of a punc- 
turing meteoroid is the surface zenith ( p n , 7^ ) . But when D < Jt/ 2 + 0, 
the probable radiant is on the earth's meridian (of the exposed hemi- 
sphere) and at ah angular increment, 7 , beyond the surface zenith, suf- 
ficiently approximated by 


xRi 

7 1 - R, 


when ■ • xRl _ g 4/3 

i. — K. -1 


= 0 + D/ 3 when 


xRi 
1 - R, 


> 4/3 


(35) 
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When D ^ re/ 2 - 0, the earth is entirely above the surface horizon, 

R 2 becomes unity in equation (31), and D replaces x in both equations (31) 
and (35). But when D is in the interval jt/2 - 0<D< it/ 2-1-0, only a 
segment of the earth* s disk is above the surface horizon, the radius of 
the earth's disk to either horizon point of the periphery makes an 
angle a with the meridian of the exposed hemisphere, and the centroid 
of the apparent segment of the earth's disk is displaced by an angle z 
from the earth* s center. Then a is a function of D and 0, z is a func- 
tion of cc and 0, R 2 in equation (31) is a function of a, and x in equa- 
tions (31) and (35) is a function of D and z. These relations are 
approximated sufficiently from the plane geometry for circular segments 
(see Reference 24) as follows: 


cos a = (D - re/ 2 ) / 0 , 0 < a < jt, (36) 

R 2 = (a - sin a cos qi)/jt, (37) 

z = (20/ 3) (s in 3 a)/ (a - sin a cos a), (38) 

x = D - z. (39) 


Let M represent the angle at the centel: of the earth's disk between 
the meridian of the exposed hemisphere and the celestial meridian; and 
let the celestial coordinates of the probable radiant of the puncturing 
meteoroid be p r and A r . Then the spherical triangles which are formed 
on the celestial sphere by the celestial pole and the great-circle arcs 
D and D + 7 have the following solution (see Reference 25): 


hav M = [hav (jt/ 2 - 6 ) - hav (jt/2 -6 - D) ] esc (tc/ 2 - 

n p 

hav (jt/2 - (3^) = hav (tt/ 2 - p^ - D - 7) 4- sin (tt/ 2 - p^) 
hav (A - A n ) cos P r cos p^ = hav 7 - hav (p R - P r ) . 


p ) esc D, 

(40) 

sin (D + 7) hav M, 

(41) 

(42) 
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Then, by equations (33), (40) , and (41), the celestial latitude, g r , of 
the probable radiant (in the interval -tc/2 ^ fa ^ jt/2) can be calculated 
from: 


sin 6 = sin(R + D + 7 ) + s , ^ n ( D 7) [ s i n « - sin(p + D)], (43) 

where, by definition, — at/ 2 ^ Pj- ^ jt/2. Also, by equations (33) and (42), 
the celestial longitude, A^., of the probable radiant can be calculated 
from 


cos (A 

\ r 


A ) = 1 + 

n 


cos y - cos (6 - 6 ) 

M n [ r 

cos 6 cos R 
! r H n 


(44) 


In addition to equation (44) for cos (A r - ’ A n ) , one must have a 
formula for sin (A r - 7^) to remove the quadrant ambiguity. Let N be 
the angle between the celestial and surface meridians at the probable 
radiant; i.e., 0 ^ N ^ it, 


hav N = [hav (jt/2 - p fl ) - hav (jt/2 - P r - 7 )] esc (jt/2 - p ) esc 7 . 

(45) 

Then, by equations (33) and (45), N can be calculated from 


cos N = 1 + 


sin R - sin (8 + 7 ) 

K n r ; 

cos sin 7 


0 ^ N £ it. 


(46) 


Therefore, by the law of sines. 


^ v sin y sin N sin 7 sin N 

iin (A r “ V - sin ( rt /2 - p") = cos p 


(47) 
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VII. CONCLUSIONS 


In appropriate meteoroid measurement satellite experiments where 
vehicle attitude information is available without information about the 
meteoroid closing velocity vector, the probable radiant and the puncture 
flux enhancement factor for sporadic meteors should be based on the 
assumption that velocity vectors for ambient meteoroids are isotropically 
distributed. Then, although half of the meteoroids which impact into a 
flat transducer are more than 45 degrees from the normal, those within 
45 degrees puncture three times more effectively than the others, and 
half of the puncturing meteoroids are less than 33.5 degrees from the 
normal when the earth is below the surface horizon. 

The probability density functions f (x 2 ) and f (x) and the correspond- 
ing cumulative distributions for angles of impact x 2 and x for incident 
and puncturing meteoroids, respectively, are as follows when the earth 
is below the surface horizon: 


f(x 2 ) = sin 2 x 2 


(18) 


/ f(x 2 ) dx 2 = (1 - cos 2 x 2 )/2 
o 


f(x) = 1.92 (cos x) 1 ’ 84 sin 2x 


(30) 


x 

J f (x) dx = 1 - (cos x) 3 * 84 . 
o 

The probability density functions in equations (18) and (30) above 
also imply that the flux 0 of puncturing meteoroids is five times greater 
than the flux F s of incident meteoroids with mass m equal to or greater 
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than the mass m of the almost-puncturing meteoroid with the average 
values of material density, velocity, and impact angle; i.e.. 


F = 10^ 6 m^ 2 
s 


( 5 ) 


0 = 


10 P6 +0 «70 


-Bp 

irr 2 . 


(29) 


The Pegasus-type experiments give puncture flux, 0 , as a function 
of target thickness, p, and target material parameter, k 2 , rather than 
as a function of the parameters p 2 an ^ Pe relating flux and mass in 
equations (5) and (29) above; i.e.. 


m = 10 kl (10 k2 p) k3 


( 12 ) 


log 0 = (p 6 + °* 70 + k lfe) + k 3fe( k 2 + l°g p) 


(14) 


where k x and k 3 depend on the unknown meteoroid material parameters, on 
the mechanics of impact, and possibly also on k 2 and p, and where k 2 is 
a known function of the known target material parameters as follows: 


k 2 = 0.50 log C t + 1.31 log E + 8.0 log + 0.43 log e^. (13) 

As observed from the site of puncture in a Pegasus- type experiment, 
the angular radius, 0, of* the earth’s effective meteoroid encounter 
cross section is sufficiently approximated from: 


sin 0 = (r g + 100)/r (16) 

where r is the geocentric radial distance of the satellite in kilometers 
and r e is the radius of the earth in kilometers. 

The separation D(0 g D £ it radians) of the center of the earth’s 
disk (at celestial coordinates Pp and Ap) from the surface zenith (at 
celestial coordinates p^ and A n ) is calculated from: 


cos D = cos (8 - 8 ) - [1 - cos (A - A )] cos 8 cos 8 . (34) 

p n n p n p 


20 



When D ^ tc/ 2 + 0 the earth is below the surface horizon and the 
probable radiant of a puncturing meteoroid is the surface zenith (Pn*^)- 
But when D is in the interval at/2 - 0<D< at/ 2-1-0, a fraction R 2 of 
the earth 1 s disk with central angle 2 oc is apparent above the surface 
horizon and shields a fraction of the meteoroids from the exposed 
celestial hemisphere. The centroid of the apparent segment is separated 
z radians from the center of the earth* s disk and x radians from the 
surface zenith. Then the probable radiant is on the earth's meridian 
(of the exposed hemisphere) and at an increment y radians beyond the 
surface zenith. The relations between these parameters are sufficiently 
approximated by: 


cos a = (D- rc/2)/0, 0 < Ctf < it 

R 2 = (a - sin GC cos a) /at 
z = (20 sin 3 a)/ (3 ttR 2 ) 
x = D - z 

R-l = 1.67(cos x ) 2,84 (1 - cos 0 ) R 2 
when - — ^4/3 

1 - R-l 

when - — >4/3 

1 " K l 


xRi 

7 " 1 - R-l ’ 


= 0 + D/3 


(36) 

(37) 

(38) 

(39) 
(31) 

(35) 


Also, when D ^ re/ 2 - 0 the earth is entirely above the surface horizon, 

R 2 becomes unity in equation (31) and D replaces x in both equations (31) 
and (35). Then, when at least part of the earth is apparent above the 
surface horizon (D < at/ 2 •+ 0), the celestial coordinates of the probable 
radiant (p r , A r ) can be calculated from: 


sin 8 = 

K r 


sin (p + D + y) 


s in (D + 2 1 r c 
sin D u 


Pn - 


s in (p 4* D) ] 




(43) 
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( 44 ) 


cos 7 - cos O - 6 ) 
cos (A - A ) = 1 + — — p 


r n 


cos 6 cos 6 
K r *n 


^ x sin y sin N 

sm (A - A ) = L — 

v r n cos 6 

K n 


where 


cos N = 1 + 


sin £ n - sin(p r + 7 ) 
cos sin 7 


0 S N g n. 


(47) 


(46) 


The formulas for the celestial coordinates of the probable radiant 
(Pr^A^-of a puncturing meteoroid have been corrected for earth-shielding. 
This correction is important near the earth even under nominal circum- 
stances; e.g . 5 when the orbit height is 500 kilometers, the radius of 
the earth* s disk (as seen from the vehicle) is 70 degrees. In that 
case, when the elevations of the earth* s center are 0 , 15, 30, 45, and 
60 degrees, the probable radiants are deflected 5, 10, 20, 39, and 80 
degrees, respectively, by the earth-shielding effect. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 
Huntsville, Alabama, November 25, 1965 
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